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The strongly-coupled phase of the quark-gluon plasma (QGP) is studied here by resorting to 
a T-matrix formulation in which the medium is seen as a non-ideal gas of quasiparticles (quarks, 
antiquarks and gluons) interacting nonpertubatively. In the temperature range under study, (1- 
5) Tc, where Tc is the temperature of deconfinement, the interactions are expected to be strong 
enough to generate bound states. The dissociation temperature of such binary bound states is thus 
computed here. The more the quasiparticles involved in the binary system are heavy, the more 
the bound state is likely to survive significantly above Tc- Then, the QGP equations of state at 
zero and small baryonic potential are computed for Nf = 2 and Nf = 2 -|- 1 by resorting to the 
Dashen, Ma and Bernstein formulation of statistical mechanics. Comparisons with current lattice 
QCD data are presented. 
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I. INTRODUCTION 


The phenomenology related to the QCD conhnement/deconhnement phase transition is 
nowadays a fascinating subject in the center of intense investigations, both experimentally 
and theoretically (see e.g. [1] for a review of the topic). During the last two decades, a new 
state of matter, the quark-gluon plasma (QGP), has been intensively studied through heavy- 
ion collisions (SPS, RHIC or LHC) and is still carefully analysed. The experiments seem to 
conclude that the QGP behaves like a perfect fluid with a low ratio viscosity over entropy 
around the critical temperature of deconhnement Tc. Therefore, this observation suggests 
that a strongly-coupled phase (called sQGP) is present in this temperature range and that 
the QCD conhnement/deconhnement phase transition is much more driven by screening 
ehects of the strong interaction. In order to correctly describe the diherent mechanisms at 
work during this phase transition, or more exactly during this crossover, a lot of theoretical 
researches (lattice QCD, phenomenological approaches...) are carried out. In particular, 
hnding the QGP equations of state (EoS) is a crucial information nowadays needed. 

The aim of this work is to investigate the thermodynamic features of the QGP by resorting 
to a phenomenological approach based on T-matrix computations. This approach has the 
advantage to allow the study of bound and scattering states of the system in a whole picture. 
Relevant results have already been established for heavy-quarkonia above [2] and also for 
glueballs in the Yang-Mills plasma [3]. Moreover, in this latter reference, the EoS of the 
Yang-Mills plasma for SU(N) and G 2 have been computed thanks to the Dashen, Ma and 
Bernstein’s formulation of statistical mechanics in terms of the S'-matrix (or T-matrix) [4]. 
Such a formulation is particularly well suited for systems whose microscopic constituents 
behave according to relativistic quantum mechanics. The QGP is indeed identified to a 
quantum gas of gluons, quarks and antiquarks, which are seen as the effective degrees of 
freedom propagating in the plasma. This assumption is actually common to all the so-called 
quasiparticle approaches [5, 6]. However, thanks to the T-matrix formulation, the strongly- 
interacting regime can also be investigated here, in which bound states are expected to still 
survive above Tc [7]. 

The paper is organized as follows. Sec. H is a summary of the approach used here and 
about which detailed explanations can be found in [3, 8]. In Sec. HI, the model parameters 
are presented and discussed. In particular, the quasiparticle bare masses are extracted from 
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the T = 0 spectrum. In Sec. IV, the binary bound state spectrum above Tc is computed 
and analysed. Then, the EoS of the deconhned phase at zero baryonic potential are studied 
for Nf = 2 and iVj = 2 + 1 in Sec. V. To hnish, an exploratory work at small baryonic 
potential is carried out in Sec. VI. All our EoS are compared to recent lattice QCD (IQCD) 
ones. Our results are hnally summarized in Sec. VII. 


II. T-MATRIX FORMALISM IN STATISTICAL PHYSICS 


A. Generalities 


The results of Dashen, Ma and Bernstein [4] establishing the grand potential of an inter¬ 
acting relativistic particle gas O, expressed as an energy density, is given by (in units where 
h = c = ks = 1) 


O — Oq + ^ ^ 




r^.E^M^E) Tr, 
/m. 47r* 


SS 




( 1 ) 


This equation is made of two parts. The hrst term Oq refers to the grand canonical potential 
of the free relativistic (quasi)particles, while the second term accounts for interactions in 
the plasma. This latter is made of a sum running on all the species, the number of particles 
included, and the quantum numbers necessary to £x a channel. The vectors jJ = (pi, /r 2 , • • •) 
and N = {Ni, N 2 ,... ) contain the chemical potentials and the particle number of each 
species taking part in a given scattering process. The set of all these channels is generically 
denoted v. As usual, the chemical potential p is the Lagrange multiplier associated with the 
number of particles. It is a measure for the density of particles. In relativistic models, the 
sign of /i is changed, passing from matter to antimatter. This parameter marks imbalance 
between matter and antimatter [4, 9]. 

One can notice that the contribution of the bound and scattering states are decoupled. 
The threshold is the summation on the masses of all the particles included in a given 
channel v. Below bound states appearing as pole in the S'-matrix (equivalently T- 
matrix) are added as free additional species: is the grand canonical potential describing 

a free relativistic gas of the //-channel bound states. Above M^, the scattering contribution 
is expressed as an integration depending on a trace, taken in the center-of-mass frame of the 
particles in the channel z/, and function of the S'-matrix of the system. S is in particular 
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a function of the total energy E. The symmetrizer S enforces the Pauli principle when a 
channel involving identical particles is considered, and the subscript c means that only the 
connected scattering diagrams are taken into account. K 2 {x) is the modihed Bessel function 

^ _ y 

of the second kind and = 1/T where T is the temperature. The symbol AdxB denotes 
A{dxB) - {dxA)B. 

By dehnition, S is linked to off-shell T-matrix T: 


^ = 1 - 27ri 6{E - Ho) T, 


( 2 ) 


where Hq is the free Hamiltonian of the system. As in [3, 8], we will only focus on two-body 
channels. So, a way to obtain T is to solve the Lippmann-Schwinger equation, schematically 
given by 

r = v + vGor, (3) 

with Go the free two-body propagator and V the interaction potential. It is worth mentioning 
that for three-body channels, Faddeev equations should be used in order to eliminate the 
spurious solution of the Lippmann-Schwinger equation [10]. Such considerations will be thus 
out of scope in this paper. 

Once Eq. (1) is computed, all thermodynamic observables can derived. For example, the 
pressure is simply given by 

p = —O. (4) 


The sum appearing in (1) explicitly reads J2jpc where only two particles are 
involved in the interaction process, / is a possible isospin channel, C is the color channel, 
and is the spin/helicity channel (the labels C or P must be dropped off if the charge 
conjugation or the parity are not dehned). 

The normalized trace anomaly can also be computed by the following formula 


A 

PSB 


-/S 



( 5 ) 


where psB is the Stefan-Boltzmann pressure. Although we give here some results about 
the trace anomaly, it is mentioned in [3] that some improvements must be done in order to 
obtain a fully reliable estimation of this quantity. 
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B. Quasiparticle properties 


Assuming that the dominant scattering processes are the two-body ones, a key ingredient 
of the present approach is the two-body potential V, encoding the interactions between 
the particles in the plasma. As in [3, 8], V is extracted from the static quenched SU(3) 
IQCD free energy Fi, between a qq pair in singlet representation at hnite temperature 

[11] , and then fitted with a Cornell potential, screened thanks to the Debye-Hiickel theory 

[12] (see Appendix B in [3]). Note that unquenched IQCD results are also available in [13]. 
Nevertheless, since these results are not signihcantly different from the quenched ones, the 
quenched potential is kept as basis of our computations, giving the accuracy expected in our 
work. 

From that, the internal energy Ui is computed, Ui = Fi — TOtFi, and considered as 
the interaction potential. This choice is still a matter of debate. Nevertheless, it has given 
correct results in the ordinary YM case, as shown in [3]. Moreover in Sec VC, we will 
see that, according to our prescription for the quasiparticle masses, the internal energy is 
required to have a better agreement between our results and IQCD ones just above Tc- 

No relativistic corrections will be taken into account for light-quark interactions within 
this paper. Indeed, the quasiparticle quark masses used in our approach are large enough 
to assume static potentials at hrst approximation. Nevertheless, this task is left for future 
works. 

Moreover, all hyperhne interactions are neglected. We can expect that they are non¬ 
dominant with respect to the spin-independent contributions, since these processes are as¬ 
sumed to depend on the inverse square of the effective mass. With this hypothesis, we also 
miss the diagonal annihilation contributions. 

Finally, the Casimir scaling is used to extract the leading-order gauge dependence of 
Ui{r,T) for T > Tc, as proposed in Sec. II in [3]. The Casimir scaling means that potentials 
between colored sources are proportional to the eigenvalues of the quadratic Casimir operator 
for their representation [14]. It is the simplest color dependence for two color sources: It 
has indeed the same form as the one for the one-gluon exchange process. Nevertheless, it is 
important to stress that the interaction considered within this paper contains other processes 
since it stems from a IQCD computation. Let us note that the annihilation mechanism, which 
does not respect the Casimir scaling, is a contact interaction and is then vanishing for all 
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non-S states. Moreover, it is worth mentioning that the Casimir scaling seems very well 
respected between two static color sources in the T = 0 sector [14], Computations in the 
T > Tc sector show a situation which is slightly different: The Casimir scaling seems partly 
violated (at most 20%) for short distances and temperatures near Tc [15]. Nevertheless in 
this work, as in [3], we assume that the Casimir scaling is satished. The hnal form of the 
potential is thus the following one: 


V(r,T) = ^[U,{r,T)-U,{oo,T)], 


K,: 




where 


^C;ij 


r*c r^Ri 

U -2 '^2 


2Cf^ 


( 6 ) 


( 7 ) 


and where C 2 is the quadratic Casimir of the representation R. C, adj, Ri and Rj stand 
respectively for the pair, adjoint, i- and j-particle representation. For instance, 

7V2 _ 1 


^ ^ Cl ^ cl = 


2N ’ 


( 8 ) 


for the SU(N) gauge group (the singlet representation is denoted by •). All the values taken 
by (7) for the various color channels considered in this study are given in Appendix A in [3]. 
Let us note that the interaction can be attractive or repulsive. The normalization of (6) is 
given by K,,.gg = —4/9, since Ui{r, T) is htted on a singlet qq potential for SU(3). We can also 
notice in (6) that the long-distance behavior of the lattice potential Ui{oo, T), is subtracted. 
Indeed, this term is assimilated, as suggested in [16], as a thermal mass contribution for 
the quasiparticles. Moreover, it ensures the convergence of the scattering equation and the 
possibility to perform the Fourier transform. 

When the quasiparticles are inhnitely separated, the only remaining potential energy 
can be seen as a manifestation of the in-medium self-energy effects, Ui{oo,T) = 2mq{T). 
We thus encode these effects as a mass shift 6(T) to the “bare” quasiparticle mass mo, by 
following the arguments exposed in [3]: 


m{TY = ml + S{Tf. 


( 9 ) 


In order to get the thermal mass for any particles, the first-order color dependence is ex¬ 
tracted in agreement with the hard-thermal-loop (HTL) leading-order behavior [17]: 





( 10 ) 
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where the quantity A(T) is assumed to be color-independent. As Ui{r,T) is htted on a 
singlet qq potential for SU{3), we have here 


^ = y^A(r) = ^A(r). (ii) 

For further details about the behavior of m{T), one can refer to Sec. V in [3] and to Sec. Ill 
in this paper. At this stage, one has to have in mind that chiral symmetry is not taken into 
account in our formalism. Comments about that issue will be given in the conclusions. 


C. Solving Lippman-Schwinger equations 


The Lippman-Schwinger equation leading to the on-shell T-matrix can be computed from 
(3) as in [3, 8]: 


Tu{E-,q,q) 


Vu{q,q) + -^j dkk‘^Vy{q,k) 

X Go{E- k) %{E; k, g') [1 ± /pi(ei)] [1 ± ^2(62)], 


( 12 ) 


where E is the energy in the center-of-mass frame, e* the asymptotic energy of the particle 
i, and where the free two-body propagator is computed thanks to the Blanckenbecler-Sugar 
(BbS) reduction scheme. Its explicit form is given in Appendix C in [8]. Moreover, the in¬ 
medium effects, namely the Bose-enhancement and the Pauli-blocking are included following 
[18]. fp is thus the distribution function of the p-species: 


/p(e) 


1 

ip 1 ’ 


(13) 


the — stands for bosons while the -|- for fermions, and p is a possible chemical potential. 
The sign choice in (12) also depends on the nature of the particles: -|- for bosons and — for 
fermions. Let us note that the impact of these in-medium effects on our EoS is very small. 
Therefore, the results obtained in [3] remain valid. 

Concerning the interaction potential IA(g, q') entering in (12), it is obtained by the Fourier 
transform of the interaction extracted in IQCD. Since our potential has a spherical symmetry, 
we have 

V{q,q,9g,g>) = 4:7rj^ dr rV (14) 
where Q = A/g^~+~g^^^-'2ggA!os"0^ and is the angle between the momenta q and q'. 
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For channels given by ordinary states, Vu{q,q') is obtained from 

yL{q,q) = 271 j dxPL{x)V{q,q,x), (15) 

where Pl is the Legendre polynomial of order L. The spin S is not indicated since our 
interaction is spin-independent. 

When at least one particle is transverse, the helicity formalism [19] has to be used. It is 
then very convenient to decompose a helicity state in the basis states in order to 

perform the computations. For a particular helicity state it reads 

L,S 

Then, it can be shown that 

V>(9.9') = Eclsn(«,9'). (17) 

L,S 

since our interaction is spin-independent. All the helicity states needed for this study are 
listed in Appendix B in [8]. 

The Haftel-Tabakin algorithm is a reliable procedure to solve the T-matrix problem 
[2, 20]. The momentum integral is discretized within an appropriate quadrature, thus turning 
the integral equation into a matrix equation, namely '^PikTkj = Vij, where schematically, 

P=l-wVG{l±fp,){l±U (18) 

and where w denotes the integration weight. The solution follows trivially by matrix in¬ 
version. Bound states are naturally poles below M^y. An interesting criterion for hnding 
them is to use the determinant of the transition function P (referred to as the Fredholm 
determinant) since it vanishes at the bound state energies [3]. Finally, once T{E-,q,q') is 
known, the on-shell T-matrix is readily obtained as T(T; q{E), q{E)), with q{E) given by 

- jmi + m 2 )^) (T 2 - (mi - 7712)^) 

2E 

III. MODEL PARAMETERS 

A. Assumptions 


Before fixing the parameters and applying the general formalism described in the previous 
section to the Nf = 2 (-1-1) QGP, let us discuss some general assumptions that we have done 



within this model. 


In our approach, there are different species of quasiparticles: the gluons {g), the light 
(anti)quarks (/, resp. 1), the strange (anti)quarks (s, resp. s) and the heavy (anti)quarks 
(c and b, resp. c and b). The (anti)quarks are spin-1/2 particles belonging to fundamental 
(conjugate) representation of the gauge group. Despite their nonvanishing mass, the gluons 
are transverse spin-1 bosons in the adjoint representation. The gluon mass is dynamically 
generated by self-energy effects, which does not imply a drastic change of their nature. It 
has been shown that the gluon must be considered as a transverse spin-1 boson to reproduce 
correctly the expected glueball spectra at T = 0 [21]. Moreover, some lattice data support 
the presence of massive transverse modes only in a gluon plasma [5]. The two-body channels 
to be considered are gg, qq, qq, qq, qg and qg. The lowest corresponding spin/helicity states 
are given in Appendix B of [8], and the possible color channels can be found in Appendix A 
of [3]. Within this study, we only focus on the SU(3) gauge group. 

Let us examine the different possibilities of interactions: 


• As explained above, the interaction (6) between two gluons or two quarks follows 
strictly the Casimir scaling and neglects all hyperhne corrections, annihilation ones 
included. 


• Although this interaction is expected to take into account complicated exchanges (since 
it stems from a IQCD calculation), it is interesting to look at the simplest possible 
Feynman diagrams between two particles. Two gluons or two quarks can exchange a 
gluon, but the basic gluon-quark interaction is a quark exchange. So the choice (7) for 
the color factor is questionable for this particular interaction. To correct this point is 
beyond the scope of this work, but it is worth mentioning that the contributions of 
the qg and qg interactions is expected to be very weak in our model (see Sec. V). 


• Processes transforming a gg pair into a qq pair exist, but we have checked that mech¬ 
anisms of order 1 are naturally suppressed since there is no overlap between gg and qq 
states [8]. As we neglect second order processes, as hyperhne interactions, we do not 
take into account transition between gg and qq pairs. 
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B. Potential at T = 0 


In order to fix our parameters for starting the computations at finite temperature, and 
to check the validity of our model, some pieces of information can be extracted from the 
T = 0 bound state spectrum as in [3]. 

In quenched SU(3) lattice QCD, the potential between a static quark-antiquark pair at 
zero temperature is compatible with the funnel form 

Vf{r) =ar-^^, (20) 

where a = 0.4 and a = 0.176 GeV^ (standard values for the running coupling constant a and 
the string tension a at T = 0). Again, we neglect the contributions of annihilation processes. 
Since the Fourier transform of Vf{r) is not dehned (because of a non-zero asymptotic value), 
a string-breaking value Vsb, has to be introduced in order to make it convergent [2]. Vsb 
is thus seen as the energy above which a light quark-antiquark pair can be created from 
the vacuum and breaks the QCD string. This scale is then subtracted and the potential 
effectively taken into account is Vf(r) — Vgb, while Kfe/2 is interpreted as an effective quark 
mass using the same arguments as those detailed in Sec. IIB. 

According to the color scaling (7), the potential describing the interactions between two 
color sources (with representations R and R) at zero temperature, is 

ro(r) = I (c« + C«) Vf(r) - Vf, (21) 


since C* = 0. The factor 9/4 appears since the potential Vf is htted on a singlet qq pair 
for a SU(3) gauge group. In this case, should rather be interpreted as the energy scale 
necessary to form two sources of color compatible with the existence of the two new color 
singlet pairs of particles created by the string breaking. As in [3, 8], if mo is the bare mass 
of the particle, the T = 0 mass m(0), used to compute the bound state is then 


m{0Y 



( 22 ) 


keeping the same structure as in (9). 
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C. T = 0 Bound State Spectrum 


The zero-temperature spectrum of the theory can be computed by solving (12) with the 
potential (21). As mentioned in Sec. IIC, instead of looking at the pole of the T-matrix, 
the zeros of det are computed in order to establish the bound state spectrum. 

The lightest glueball spectrum, namely the O’*"’*', 0 and 2+’'', has already been computed 
in [3]. The parameters, and mo, were respectively hxed to 2 GeV and 0.7 GeV. 
is in agreement with lattice data showing that the mass of the lightest gluelump is given 
by 0.85(17) GeV [22], while mo is an acceptable value for the zero-momentum limit of the 
gluon propagator at zero temperature in view of previous studies locating this mass typically 
between 500 and 700 MeV (see e.g. [23-25]). Interesting reader can refer to [3] for additional 
information about the T = 0 glueball spectrum. 

Within this paper, the stress is put on mesons with an orbital angular momentum L = 0 
or L = 1. The allowed states with these quantum numbers are displayed in Table I. As it 
can be seen, several states are associated to a same L. Since the potential (21) does not 
depend on other quantum numbers, all these states are degenerate within our approach. 


J 

L 

S 


0 

0 

0 

0-+ 


1 

1 

0++ 

1 

1 

0 

1+- 


0 

1 

f— 


1 

1 

1++ 

2 

1 

1 

2++ 


TABLE I. states allowed for qq with L = 0 or L = 1 at T = 0. The parity of the state P, is 
given by (—l)^"*"^ while the charge conjugaison C, is (—l)^"’"'®'. 


The used parameters are summarized in Table II. There are essentially two main points to 
notice. Firstly, a shift of 0.3 GeV to the PDG quark bare mass [26] is systematically present. 
It is a common assumption within quasiparticle approaches since this shift corresponds to 
one third of the nucleon mass. Moreover, it is a typical value for the chiral condensate 
according to [27]. Secondly, the string breaking depends on the quark flavor. This could 
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be explained by the following argument: It is not the same region of the potential that is 
relevant for the dynamics of all the quark flavors. Indeed, light quarks are more sensitive 
to the linear part of the interaction while the heavy-quark potential is dominated by the 
Coulomb one. 

According to this interpretation, has to be higher for light quarks and has to decrease 
with heavier quark flavors. It is exactly what it is observed with the parameters in Table II. 
Moreover, when the bound state is made of two different quark flavors, takes a value 
between the chosen string breaking for the two associated quarkonia systems. For D and 
5-mesons, is closer to the string breaking of heavy quarkonia. 


Quark composition 

vS 

ml 

ml 

Light (1-1) 

2.6 

0.3 

0.3 

Strange (s-s) 

2 

0.4 

0.4 

Charm (c-c) 

1 

1.6 

1.6 

Beauty (b-b) 

0.7 

4.95 

4.95 

Kaon (1-s) 

2.4 

0.3 

0.4 

D-meson (1-c) 

1.5 

0.3 

1.6 

Ds-meson (s-c) 

1.2 

0.4 

1.6 

B-meson (1-b) 

1.2 

0.3 

4.95 

B^-meson (s-b) 

1 

0.4 

4.95 

Bc-meson (c-b) 

0.7 

1.6 

4.95 


TABLE II. Masses and string breaking (in GeV) for the different flavors of quarks. 

In Table III, the results are compared to experimental [26]. As it can be noticed, a 
quite good agreement is reached provided that we do not consider the lightest pions and 
kaons i.e. 7r(140) and iF(495). Indeed, the fact that the mass of these lightest mesons are 
not achievable can be explained by the theoretical origin of such states: The pion is the 
Goldstone boson resulting from the spontaneously breaking of the chiral symmetry. A so 
peculiar phenomenon can not be described within such simple effective model. Moreover, 
according to quasiparticle standard approaches, the spin effects are the weakest in a 5 = 1 
channel. Since our computations do not take into account such effects, it is reasonable that 
our results for the L = 0 light mesons are closed to the p instead of the tt. 
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L = 0 

Exp. 

T -matrix 

L = 1 

Exp. T-matrix 

p{uu, dd) 

0.77 

0.72 

ao{uu, dd) 

1.45 

1.45 

$(ss) 

1.02 

1.08 

f2iss) 

1.53 

1.58 

K*{l-s) 

0.89 

0.89 

K{l-s) 

1.43 

1.52 

d/V'(cc) 

3.10 

3.01 




r{bb) 

9.46 

9.40 




D*{l-c) 

2.01 

2.01 




Dl{s-c) 

2.11 

2.11 




B*{l-b) 

5.33 

5.33 




B*s{s-h) 

5.42 

5.41 




B*{s-c) 

- 

6.39 





TABLE III. Masses (in GeV) of the L = 0 and L = 1 meson states at zero temperature with the 
gauge group SU(3). Our results (third and sixth column), are compared to the experimental data 
of [26] (second and fifth column). 


Finally, let us add that, unlike in the glueball case, the T = 0 meson mass depends on the 
gauge group since K,,.gg depends on it (see Appendix A in [3]). Within our approach, such 
study is not difficult to carry out. The meson mass dependence in function of the gauge 
group is not studied here, since the principal interest of the T-matrix computations at T = 0 
is to extract and to check the parameters we will use at T 7 ^ 0. In this regard, let us note 
that only the quasiparticle bare masses will enter in our computations at T 7 ^ 0 and not the 
string breaking Vgb- What can be said however about the gauge-group dependence is that, 
in the case of SU{N), the meson masses are of order 1 as expected, see conclusions for more 
comments about this. 


D. Critical temperature of deconfinement 

The last global parameter that has to be hxed within our approach is the critical temper¬ 
ature of deconhnement Tc- As we £t our interactions on lattice calculations, the definition 
of Tc comes from these approaches: the color averaged as well as the singlet free energy of 
a quark-antiquark pair will tend towards finite nonzero values for all temperatures T > 
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(and diverge below Tc). In [3], a value of 0.3 GeV was used since the focus was only on 
the gluon sector. Here, in order to stay coherent with current lattice data, is moved to 
0.15 GeV. This change naturally modihes the thermodynamics in the gluon sector that was 
established in [3]. Let us discuss this point. 

First of all, let us set = T/Tc. The two-body interaction potential between particles 
only depends on as it can be explicitly shown from its expression given in Appendix B of 
[3]. Therefore, it is the same for the gluon thermal mass Sg{T), according to (10). Moreover, 
since the effect of the Bose-enhancement can be considered as negligible, it can also be 
assumed that the T-matrix Tjpc, is only a function of 2 :. 

In Fig. 1, the Tc-impact is analyzed for the pure gauge EoS with Tc = 0.15 and 0.30 GeV, 
all other parameters remaining hxed. The way to compute these EoS is given in [3] and will 
be recalled in Sec. V A. As it can be noticed, the general behavior is not the same and the 
normalized pressure seems to increase when Tc increases. It is especially worth remarking 
that, if we only consider the contribution of the quasiparticle ideal gas, Tc has an impact 
on the thermodynamics. Indeed, the pressure depends on the ratio m/T = m{z)/{zTc), 
depending on Tc for a fixed 2 :. 



FIG. 1. Normalized pure-gauge pressure p/psB, and trace anomaly A/pss (without bound states) 
versus temperature in units of Tc (with Tc = 0.15 and 0.3 GeV). 


The behavior of the normalized trace anomaly (without bound states) is also presented 
in Fig. 1. The Tc-dependence is not very easy to predict because the trace anomaly depends 
on the slope of the associated pressure curve, and small variations can generate a drastic 
change. In Fig. 1, we can indeed observe that the behavior around Tc is extremely different 
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for Tc = 0.15 and 0.3 GeV. The peak structure is lost for Tc = 0.15 GeV. This could be 
due to two reasons: the total change of the free part structure and the small impact of the 
interactions in comparison with the EoS obtained in [3], as it is shown in Fig. 2. 



TIT, 



FIG. 2. (Left) Normalized pure-gauge pressure p/psB versus temperature in units of Tc (with 
Tc = 0.15 GeV), compared to the free part contribution. (Right) Normalized pure-gauge trace 
anomaly A/pss versus temperature in units of Tc (with Tc = 0.15 GeV), compared to the free part 
contribution. 


IV. BOUND STATES WITHIN THE QGP 

The existence or not of bound states in the deconhned phase is not forbidden in principle, 
especially around Tc where interactions are expected strong enough to bind two or more par¬ 
ticles [7]. Since the operator Kc-ij is negative for several color channels, the hnite-temperature 
spectrum of QGD above Tc can be computed by solving (12) with the potential (6). The 
thermal masses of the quasiparticles are given by (9) with tuq extracted from the T = 0 
bound state spectrum (see Sec. IIIG). 

Within our formalism, the channels in which bound states are favored at most should 
contain a S'-wave component to avoid the centrifugal barrier and should have a symmetry 
that allows the state to be in a color singlet, the color channel in which the interactions are 
maximally attractive. 

In the gg case there are two such states: The O’*"’*' and 2~^~^ ones, in color singlet, cor¬ 
respond to the scalar and tensor glueballs respectively. We have observed in [3] that both 
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Light quark sector 

Light-strange quark sector 

Channel 

qq {L = 0) 

qq/qq {L = 0) 

Channel 

qq {L = 0) 

qq/qq (L = 0) 

C 

Singlet 

AS 


C 

Singlet 

AS 

T/Tc 

2 mi 



T/Tc 

mi + ms 



1.05 

1.67 

1.51 

1.67 

1.05 

1.71 

1.54 

1.71 

1.10 

1.28 

- 

- 

1.10 

1.43 

- 

- 


Strange quark sector 


Channel 

C 

qq {L = 0) 

Singlet 

qq/qq {L = 0) 

AS 

T/Tc 

2 ms 



1.05 

1.76 

1.57 

1.74 

1.10 

1.38 

1.38 

- 

1.15 

1.23 

- 



TABLE IV. Masses (GeV) of lowest-lying QCD spectrum above Tc {Tc = 0.15 GeV). Singlet and 
AS respectively refer to the singlet and antisymmetric representation of the gauge group. A line 
mark the temperature at which a bound state is not detected anymore. 

the scalar and tensor glueball masses at 1.05 were compatible with the zero-temperature 
ones. Moreover, the scalar glueball exists as a bound state up to 1.25 while the tensor 
one is bound up to 1.15 Tc- Note that in [3], the impact of the Bose-enhancement were 
not considered in the T-matrix. Therefore, the value of Tc should modify the masses of the 
bound states. Nevertheless, it has been numerically checked that the data shown in [3] differ 
from the ones containing in-medium effects only with a relative error of the order of 2%. 
This is the reason why they are not presented again here. 

Concerning the light and strange mesonic sector, the dissolution inside the plasma is 
much more rapid. Indeed, the most attractive channel, the L = 0 in singlet, is the only one 
that survives above T^. Moreover, mesons quickly dissolve inside the plasma as it can be 
observed in Table VI, and the meson masses are not compatible with the T = 0 one as in 
the gg case. 

The same assertion can also be drawn for the qq and qq sector. Indeed, since there is 
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no more confinement (i.e. only singlet representation) in the QGP, such states could exist 
above T^; the most attractive one being the L = 0 in the antisymmetric representation (AS). 
Nevertheless, they also rapidly disappear just above Tc, as shown in Table VI. Such rapid 
dissolution in comparison with the gg case can be understood by the fact that the quark 
quasiparticle mass is lighter than the gluon one and that the K,-gg is more than two times 
the maximum magnitude of kc in the quark sector. Indeed, n,-gg = —1 while K,-qq = —4/9 
and KAs-,qq = i^ASm = —2/9 for a SU(3) gauge group (see Appendix A of [3]). Note that 
similar comments can also be done about the qg and qg sector, leading to a quick melting 
of these bound states inside the plasma (see Table IX). 


Charm quark sector 

Channel 

qq {L = 0) 

qq/qq {L = 0) 

qq {L = 1) 

qq/qq (L = 1) 

C 

Singlet 

AS 

Singlet 

AS 

T/Tc 

2 nic 





1.05 

3.56 

3.14 

3.44 

3.48 

- 

1.10 

3.39 

3.20 

3.36 

- 


1.15 

3.33 

3.22 

- 



1.25 

3.28 

3.26 




1.35 

3.26 

- 





TABLE V. Masses (GeV) of lowest-lying charmonium states above Tc {Tc = 0.15 GeV). Singlet 
and AS respectively refer to the singlet and antisymmetric representation of the gauge group. A 
line mark the temperature at which a bound state is not detected anymore. 

Concerning the heavy quark sector, quarkonia have already been studied within a similar 
T-matrix approach as the one proposed here [2]. The main differences are the inclusion of a 
relativistic correction to the potential in [2] and the way of implementing the quasiparticle 
masses. Within this paper, the procedure (9) to determine the quasiparticle masses is 
applied and allows one to compute systematically a large panel of binary bound states made 
of different quark flavors. 

Our study for the heavy quarkonia is displayed in Tables V and VI. Around Tc, the J/ijj 
and T masses are compatible with the T = 0 spectrum, unlike for the light and strange 
mesons. Moreover, they signihcantly survive above Tc, even if the dissociation temperatures 
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Beauty quark sector 

Ghannel 

qq (L = 0) 

qqjqq {L = 0) 

qq {L = 1) 

qq/qq {L = 1) 

C 

Singlet 

AS 

Singlet 

AS 

T/Tc 

2 rufe 





1.05 

10.2 

9.35 

9.75 

9.60 

9.90 

1.15 

9.94 

9.62 

9.87 

9.79 

9.94 

1.20 

9.93 

9.71 

9.86 

9.92 

- 

1.25 

9.93 

9.70 

9.90 

- 


1.30 

9.92 

9.75 

9.90 



1.50 

9.92 

9.79 

- 



2.00 

9.91 

9.88 




2.40 

9.91 

- 





TABLE VI. Masses (GeV) of lowest-lying bottonium states above Tc {Tc = 0.15 GeV). Singlet and 
AS respectively refer to the singlet and antisymmetric representation of the gauge group. A line 
mark the temperature at which a bound state is not detected anymore. 


are lower than the ones found in [2]. qq and qq states can also be formed with the medium 
but they dissolve more rapidly than the associated quarkonia, due their weaker interaction 
potential. 

An analyse for the D and B mesons has also been carried out as well as for the qg and 
qg states for all the quark flavor considered here. The different temperatures of dissociation 
are displayed in Tables VII, VIII and IX. We can notice that the more the quasiparticles 
considered in the binary state are heavy, the more it survives significantly above T^. 

In [28], it is found that radially excited states are unlikely to survive above Tc since 
they tend to melt below the phase transition because of string breaking effects at hnite 
temperature. In our case, such states can be found but not in all channels. However, they 
quickly dissolve within the medium. 
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Dissociation temperature: qq 


Light 

Strange 

Charm 

Beauty 

Light 

1.10 ± 0.05 

1.10 ± 0.05 

1.15 ± 0.05 

1.20 ± 0.05 

Strange 


1.15 ± 0.05 

1.15 ± 0.05 

1.20 ± 0.05 

Charm 



1.35 ± 0.05 

1.60 ± 0.1 

Beauty 




2.4 ± 0.1 


TABLE VII. Temperature of dissociation in units of Tc for L = 0 mesons. 


Dissociation temperature: qq and qq 


Light 

Strange 

Charm 

Beauty 

Light 

1.10 ± 0.05 

1.10 ± 0.05 

1.10 ± 0.05 

1.10 ± 0.05 

Strange 


1.10 ± 0.05 

1.10 ± 0.05 

1.10 ± 0.05 

Charm 



1.15 ± 0.05 

1.25 ± 0.05 

Beauty 




1.50 ± 0.05 


TABLE VIII. Temperature of dissociation in units of Tc for L = 0 qq and qq states. 


Dissociation temperature: qg and qg 

Light 

Strange 

Charm 

Beauty 

1.10 ± 0.05 

1.10 ± 0.05 

1.20 ± 0.05 

1.25 ± 0.05 


TABLE IX. Temperature of dissociation in units of Tc for qg and qg states . 

V. EQUATION OF STATE OF THE QGP AT /x = 0 
A. General expression 

Now that the bound-state sector is analyzed, it is possible to compute explicitly the EoS 
and so, to study the QGP thermodynamics. In what follows, the heavy quark states will be 
not included in our EoS. Indeed, their contributions to the bound-state and scattering parts 
of the grand canonical potential are expected to be small because of their large bound-state 
masses. 
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Some preliminary lattice results about the influence of charm quarks on the EoS can 
be found in [29]. It appears actually that charm quarks bring a significant contribution 
to the trace anomaly above T^. However, a technical problem of the present approach is 
that discontinuity appear in the trace anomaly when bound states melt [3]. This problem 
is especially apparent when heavy quarks are involved. Hence, including heavy flavors in 
our computations would lead to results that are probably not reliable, and we prefer not to 
consider them. 

Let us thus particularize (1) to a QGP with Nf = 2 (+1). As in [3, 8], a two-body 
restriction is used: The considered interactions are gg, qq, qq, qq, qg and qg, in different 
colour and channels. Therefore, the first term in (1), i.e. the free relativistic gas is 

given by 

Nf Nf 

= 2 dim adj (m^ , 0)^ + 2 ^ dim Uq {mq„ , ) + 2 ^ dim uJq {mq„ , -/ig„ ) , 

, n=l n=l 

gluons ^ ^ ^ ^ ^ ^ 

quarks antiquarks 

(23) 

where the gluons have a mass nig, and the (anti)quarks (g„) a mass given by the 
prescription (9) with the mo value extracted from the T = 0 spectrum. is the chemical 
potential of the considered quark flavor. They are set to zero within this section. So, 
no asymmetry between quarks and antiquarks is taken into account within the QGP. The 
particle degrees of freedom are the following. The gluon is a transverse spin-1 (so, two spin 
projections) boson lying in the adjoint representation of the gauge group, while the quark 
(resp. antiquark), existing in Nf different flavors, is a spin-1/2 fermion belonging in the 
fundamental (resp. conjugate) gauge-group representation. The grand canonical potential 
per degree of freedom associated to a bosonic species and to a fermionic species 

(jJq {m, n), with mass m are given by 

1 / \ 

^ ^ i, * , (24) 

^^■(m, #1) = -^ * e hi (l + e-4(vWJ -,.)y (25) 

For later convenience, the thermodynamic quantities will be normalized to the Stefan- 
Boltzmann pressure, which is defined as 

PSB = - lim (26) 

m^O 
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and reads in this case 


PSB 


71 


45/3^ 


Nf 


dim adj + - dim t 


n=l 


(27) 


As already mentioned, the second term of (1) stands for the interactions. The sum 
now explicitly reads J2ng+nq +nq =2 Xlc 'TlijP) where riq^ are respectively the number 

of gluons, quarks and antiquarks involved in the interaction process. Attractive interactions 
can lead to the formation of bound states with masses M^jp < mi +m 2 (see Sec. IV). They 
contribute also to the grand potential as new species via the formula 


ng+nq^+nq^=2 JP C 


All and color channels leading to bound states are included in this summation. 

Concerning the scattering term, a tedious calculation (explained in [3, 8]), leads to the 
following result. 


^QCD^ 


647r®/d^ 


E 


3/3(mi+/22) 


nq+nq„+nq„=2 


^(2J + 1) ^ dimC 
jp c 


(29) 


/d / dee^ uj{e) A{e)Ki[lde)'ReTc jp{e,u{e),u{e)) 


J m\+m2 

1 r 


1 


'mi+m 2 

poo 


de uj{ef k{efK2{l3e) ReTc,jp(e; a;(e), a;(e)) (imTc^jp (e; a;(e), a;(e)))' 
de u{efk{efK2{(3e) (ReTc,jp(e; a;(e), a;(e)))' ImTc,jp(e; a;(e), a;(e)) 


'mi+m2 

where u{e) and A(e) are given by 


u{e) 

A{e) 


a/( e^ - (mi + m 2 )^)(e 2 - (mi - 7712 )^) 

Ye 

— {ml — ml) 

^3 ’ 


(30) 

(31) 


and where Tcjp{e] u{e),u{e)) is the on-shell jp-matrix. Note that a isospin number I has 
to be taken into account when one deals with u and d quarks since they have the same mass 
in our approach. This isospin number enters in the summation +n-=2 ^ 

factor. 

Finally, the grand canonical potential (reduced to two-body interactions) is summarized 
by the following formula 
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figCD _ ^QCD ^ ^QCD ^ ^QCD_ (32) 

For obvious numerical reasons, the summation over the number of particles is not the only 
one that must be restricted. All possible color channels are included, but all the possible 
jP(c) channels contributing to can not be included since their number is inhnite. So, 

a reliable criteria to select the most signihcant ones has to be established. The basic idea, 
already proposed in [8], is that only states with low L are included since they are the most 
likely to contribute signihcantly to a total mean cross section ctjp. Are only retained, the 
channels for which the value of ajp is at least 25% of the value djp for the channel with the 
lowest value of (see Appendix D of [8] for further explanations). 

In the present case, this criterion implies that only the following channels are 

included. 

• For gg channels: the 0 2^^^ and I’*"’*' states; 

• For gg, gg and gg channels: all the ones with L = 0 or L = 1; 

• For the qg and qg channels: all the states with (L^) < 8 (see Appendix B in [8]). 

B. QGP with Nf = 2 

Now that the number of channels for each two-body interactions is hxed, the EoS can 
be computed. In Fig. 3, the normalized pressure is shown for a QGP with two light-quark 
flavor included. As it can be noticed, interactions do not practically contribute: The major 
part is given by the free gas. Globally, the weakness of interactions can be interpreted in the 
same way as what it is observed for the YM plasma in Sec. Ill D. Indeed, when the critical 
temperature decreases, the interactions seem to become smaller and smaller. This behaviour 
is driven by the Bessel functions entering in the dehnition of the scattering part. Moreover, it 
is important to notice that the integration range in this term formally starts at /3{mi + 7712 ). 
This value is large in comparison to the values at which the Bessel functions is signihcantly 
non-zero. A change of the thermal mass prescription could thus impact the contributions of 
the scattering part. This work is left for further developments of our approach. 

In the right panel of Fig. 3, the different scattering contributions are separated. Without 
surprise, the gg, qg and qg channels asymptotically tend to zero. Indeed, it has been shown in 
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FIG. 3. (Left) Normalized pressure p/psB versus temperature in units of Tc, compared to the free 
part, bound state and scattering contribution. (Right) Different scattering contributions to the 
normalized pressure p/psB versus temperature in units of T^. Tc = 0.15 GeV. 


[3] that the interactions between two different species vanish within the Born approximation, 
because of an identity relating the color factors: Xlc dimC Kc,ij = 0. Concerning the qq and 
qq channels, they generate a global increase of the normalized pressure while it is the contrary 
for the gg sector. Not only these two effects are weak but in addition, they contribute in 
opposite directions, leading to a global suppression of the two-body interactions in average. 



TIT, 


FIG. 4. Normalized trace anomaly A/psB (without bound states) versus temperature in units of 
Tc, compared to the free part contribution. Tc = 0.15 GeV. 

In Fig. 4, we display the normalized trace anomaly (without bound states) compared 
to the free gas part. A peak structure is here exhibited even in the free gas contribution. 
Therefore, it is different from the YM sector where the interactions create the peak. The 
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nature of this latter is really difficult to establish since few variations of the pressure can 
drastically change the shape of the trace anomaly. 

The main conclusion that seems to emerge from our approach (looking at the normalized 
trace anomaly as well as at the normalized pressure) is that the leading behavior of the QGP 
is driven by gluon and (anti)quark degrees of freedom that interact weakly. Nevertheless, it 
does not mean that the interactions have no impact on the EoS. Indeed, the particle thermal 
mass is extracted from it, leading to a self-energy contribution for the particle (see Sec. V C). 


C. QGP with Nf = 2 


A similar analyse as the one proposed in the previous subsection can be applied in the 2 
-|- 1 QGP case. Since similar results and features can be deduced from it, we will not repeat 
it again and focus more on the comparisons between our model and IQGD extracted from 
[30]. Indeed, IQGD collaborations have recently reached the physical quark masses in their 
computations of the EoS, making their results more and more reliable for comparisons. 




FIG. 5. (Left) Normalized pressure p/psB versus temperature in units of Tc, compared to IQCD 
data from [30]. (Right) Normalized trace anomaly A/psB versus temperature in units of Tc, 
compared to IQCD data from [30]. Tc = 0.15 GeV. 


As we can observe in Fig. 5, our data are qualitatively in agreement with IQGD ones. The 
IQGD normalized pressure is slightly overestimated as well as the asymptotic behavior of 
the normalized trace anomaly. On the other hand, the peak structure of the trace anomaly 
is very different of the IQGD one. As already mentioned, this latter is really difficult to 
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obtain due to several reasons in our approach: problems in the inclusion of the bound state 
[3], reliability of the quark masses and restriction to two-body interactions. Even in IQCD, 
different collaborations hud different quantitative behaviors for the trace anomaly peak [30]. 
The possible discrepancies can arise from the choice of the fermionic lattice action, the 
lattice spacing, the considered quark masses, the extrapolation to the continuum limit... 
Only, a good agreement in the behavior of the decreasing tail is reached by the different 
IQCD groups, according to [30]. Nevertheless, it is worth mentioning that the disagreements 
observed in the quantitative value of the peak structure in various IQCD results (mainly due 
to a computation with no physical quark masses) seem to reduce, and the shape of the IQCD 
trace anomaly tends to the one depicted in Fig. 5, and hrstly given by the BMW collaboration 
[31]. 

Therefore, except for the normalized trace anomaly peak structure (for which a more ap¬ 
propriate treatment of the bound-state inclusion is needed and could change signihcantly its 
structure), our data are in correct agreement with IQCD ones. As discussed in the previous 
subsection, this agreement seems to be reached by only including a quasiparticle thermal 
mass: The contributions of the two-body interacting channels are minor. Nevertheless, it is 
worth insisting on the fact that the thermal mass effects are extracted from the two-body 
IQCD interaction potential within our model. So, the chosen two-body interactions are not 
useless to understand the behavior of the QGP around Tc. Indeed, if we change the po¬ 
tential, the free gas contribution is modihed since the quasiparticle thermal masses depend 
on it. This leads to a completely different behavior of the EoS as seen in Fig. 6 in which 
the potential is now chosen to be the free energy. We can especially notice in Fig. 6 that a 
better agreement between our model and IQCD normalized pressure is reached around 
thanks to the internal energy while the discrepancy between the two curves decrease when 
the temperature increases. 

Finally, let us compare in Fig 7 the normalized pressure and trace anomaly for a QGP 
with Nf = 2 and Nf = 2-|- 1 (T^ = 0.15 GeV) to the ones of the YM plasma (T^ = 0.3 GeV). 
We can notice that the normalized pressure curves are almost superimposed and that the 
decreasing trend of the trace anomaly is nearly the same in all the considered theories. The 
maximum of the deviation between these curves is around 1.2 Tc, at the localisation of the 
trace anomaly peak. It is nevertheless important to remember that the critical temperature 
and the normalization are not the same in all the EoS (see (27)). However, within these 
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FIG. 6. (Left) Normalized pressure p/psB versus temperature in units of Tc- (Right) Normalized 
trace anomaly A/pss (without bound states) versus temperature in units of Tc. In the two figures, 
the gray (black) line is the free part contribution of computed with the free (internal) energy. 

Tc = 0.15 GeV. 


units, a universality at large temperature (> 3Tc) seems to emerge. 




FIG. 7. (Left) Normalized pressure p/psB versus temperature in units of Tc- (Right) Normalized 
trace anomaly A/psB (without bound states) versus temperature in units of Tc- Tc = 0.15 GeV 
for QGP and Tc = 0.3 GeV for YM. 


VI. EQUATION OF STATE OF THE QGP AT SMALL p 

Now that the EoS for the QGP are computed and favourably compared with IQCD, 
we can investigate the non-zero baryonic regime. This latter deserves a lot of interests, 
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especially in the area of the neutron star physics. Indeed, since pioneering works [32] about 
the existence of a deconhned phase in QCD, it was assumed that the core of the heaviest 
neutron star should be probably hlled by a medium with a high nuclear density and in which 
the signihcant degrees of freedom should be the quarks. Therefore, getting the QCD EoS at 
hnite /i could shed some light in this held. 

Up to now, this task still remain difficult from hrst QCD principles. Remember that even 
in IQCD some conceptual troubles appear (cfr. sign problem) and only perturbations around 
/r = 0 are meaningful. Therefore, it seems appealing to check whether or not quasiparticle 
approaches could help. Unfortunately at the present stage, some problems also appear in 
our formalism. The main reasons are the following. 

First, the Dashen Ma and Bernstein formalism that we have used to compute the EoS is 
based on a virial expansion in terms of . We are thus limited by construction to small 
baryonic potential. Indeed, increasing the baryonic potential is the same as increasing the 
density of particles: The many-body interactions are more and more likely to contribute. 
So, the reduction to two-body interactions becomes a poor approximation a priori and some 
problems, other than a careful computation of all the channels, arise. Let us mention for 
instance, the absence of a helicity formalism for many-body systems in a potential approach, 
and the necessity to resort to Faddeev and higher equations for more than two-body interac¬ 
tions. Moreover, when the density of particles increases, the notion of quasiparticle becomes 
more and more questionable. 

Another peculiar problem is the building of a coherent interaction in presence of baryonic 
potential. Already at two bodies, no IQCD data are available to our knowledge. It is not 
only important to define the potential between particles but also the quasiparticle mass, 
which seems to rule the main behavior of the EoS at /i = 0. A way to circumvent this 
problem could be to use the HTL expressions for the particle thermal mass, but it was not 
the bias adopted within this study. Indeed, the actual shape of our thermal masses are not 
the ones extracted from HTL. 

For all these reasons, the study that follows will be only limited to small baryonic po¬ 
tentials. We will thus keep the restriction to two-body interactions which can make sense 
in such a /r-range. Moreover, the interaction potential and the quasiparticle thermal masses 
are the same as the ones used up to now, without the inclusion of the baryonic potential. 
Of course, the obtained results must be considered as preliminary and are just intended to 
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draw a general tendency. The baryonic potential enters at two levels in our computations: 
in the T-matrices because of the in-medium effects, and in all the EoS contributions as 
multiplicative factors. Fortunately as for the Tc-impact, it seems that the /i-dependence on 
the T-matrix calculations is negligible (see Fig. 8). Therefore, these latter do not have to 
be recomputed at each fx, which drastically reduces the computational time. 



FIG. 8. Real part of the T-matrix {qq in T = 0, here given as example) in function of the energy 
for T = 0.105 GeV (left), T = 0.150 GeV (middle) and T = 0.300 GeV (right) at different /U (MeV) 
with Tc = 0.15 GeV. 


In Fig. 9, we have plotted the normalized pressure and trace anomaly at different /r for 
a QGP with two light quarks. The normalization is given by (27), that is to say at m = 0 
and yU = 0. Naturally, the gluon chemical potential is zero and the quark one is such that 
Tn = Td = T- We can see in this figure that the normalized pressure increases with /r. 
This pressure is especially driven by the increase of the free quark gas contribution given in 
Fig. 10. Indeed, as in the /r = 0 case, the leading contributions to the normalized pressure 
are the free part ones since the impact of the interactions is small as observed in Fig. 11. 
Moreover, the decrease of the free antiquark gas contribution is slower than the increase of 
the free quark gas one, explaining the total increasing behavior of the normalized pressure. 

Concerning the normalized trace anomaly, it is much more difficult to understand the 
/i-dependence. The only assertion that we can do is that the trace anomaly peak becomes 
higher and higher with the increase of /i. Moreover, we can notice that the convergence to 
zero is faster with large /i. 

As already mentioned, the scattering contributions are small. Nevertheless in Fig. 11, we 
can observe a significant dependence in terms of /i. The qq and qg scattering contributions 
obviously increase with /i, respectively as and while the qq and qg ones go in opposite 
way. However, as in the free gas case, the increase is higher than the decrease, leading in fine 
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FIG. 9. (Left) Normalized pressure p/psB versus temperature in units of Tc at different p (MeV) 
for a QGP with two light quarks. (Right) Normalized trace anomaly A/pss (without bound 
states) versus temperature in units of Tc at different p (MeV) for a QGP with two light quarks. 
Tc = 0.15 GeV. 




FIG. 10. (Left) Free quark gaz contribution to the total normalized pressure p/psB versus tempera¬ 
ture in units of Tc at different p (MeV) for a QGP with two light quarks. (Right) Free antiquark gaz 
contribution to the total normalized pressure p/psB versus temperature in units of Tc at different 
p (MeV) for a QGP with two light quarks. Tc = 0.15 GeV. 


to a more important contribution of the scattering parts to the total normalized pressure. 
To be complete, the qq scattering contribution is stable since there is no /r-dependence at 
the level of the EoS: Indeed, we have = 1 and just a very weak dependence on p 

appears in the T-matrix. Moreover, since the gg sector is independent of /i, increasing p 
means increasing the impact of the quark sector within the QGP. 

Finally, we close this study by comparing our preliminary results to the IQCD ones 
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FIG. 11. (Left, top) (/^-scattering contribution to the total normalized pressure p/psB- (Right, top) 
Same for qq. (Left, middle) Same for qg. (Right, middle) Same for qg. (Bottom) Same for qq. All 
the scattering contributions are presented versus temperature in units of Tc, with Tc = 0.15 GeV, 
at different p (MeV) for a QGP with two light quarks. 


given by [33]. Within this paper, they deal with a QGP with Nf = 2 + 1 and with a 
small baryonic potential ps- Each flavor of quarks is considered to carry one third of ps- 
Therefore, we analyse the 2 + 1 QGP with Pu = IJ^d = As for the + = 0 case, we 

sightly overestimate the normalized pressure and we miss the peak of the normalized trace 
anomaly. Therefore, it seems that these differences have mainly the same origin as at /i = 0 
and our extrapolations at small p is compatible with IQGD. 
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FIG. 12. (Color online) (Left) Normalized pressure p/psB compared to lattice QCD [33] versus 
temperature in units of at different pB (MeV) for a QGP with Nf = 2 + 1. (Right) Normalized 
trace anomaly A/psB compared to lattice QCD [33] versus temperature in units of Tc at different 
Pb (MeV) for a QGP with Nf = 2 + 1. Tc = 0.15 GeV. 
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VII. CONCLUSIONS 


The present work is part of a program aiming at studying the thermodynamic properties 
of gauge theories in the deconhned phase. The interested reader may read [3] and [8] for 
pure YM and SUSY YM theories, while this paper is devoted to the “realistic” quark- 
gluon plasma. The framework developed is based on a T-matrix formulation of statistical 
mechanics, in which the thermal masses and two-body interactions are derived from the 
static potential between fundamental color source computed in quenched lattice QCD [11]. 
Apart from the potential, the only remaining parameters are the value of Tc and the bare 
quark masses. These masses are fitted on the meson spectrum at zero temperature. The 
main assumption underlying our model is actually that a quasiparticle picture of deconfined 
matter just above deconhnement is relevant. Although it is not a rigorous proof, the nice 
agreement between our computed equations of state an the recent lattice data of Refs. 
[30, 33] can be seen as an a posteriori validation of our framework. 

We are now in position of summarizing some of the key results obtained in this paper: 

• Both the free energy or the internal energy could be used as potential terms in our 
model. It appears that, keeping the same procedure and the same values for the 
parameters, only the internal energy is able to generate an equation of state which 
has the qualitative features of the lattice equation of state. The internal energy thus 
appears as the most relevant potential within in our framework and leads to a good 
agreement with the lattice equation of state. Note that this problem is far from being 
elucidated, see for example the recent work [28], where the opposite conclusion is 
reached. 

• Between 1 and 2 Tc, color interactions are strong enough to create mesons, i.e. a 
quark-antiquark bound state in a color singlet. Mesons made of one or two light 
quarks are almost all dissociated in T^. Only mesons made of two heavy quarks (c, b) 
are bound enough to survive in the range (1.3 — 2) Tc- Although we use a T-matrix 
formulation as well, our parameters have different values of that used in [2] where the 
main goal was to reproduce mesonic correlators computed on the lattice and not the 
equation of state. In this last work, the J/ip meson is bound up to 3 Tc and the T 
meson is bound up to 3.5 Tc, thus at much higher T than what we hnd. It is worth 
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recalling that we are able to compute T-matrices in channels where the quark and the 
antiquark have different masses, which was not considered in [2], 


• Although strong in the color singlet channel, the contribution of two-body interactions 
to the equation of state is weak with respect to the free-gas part. This is partly due 
to a cancellation between attractive and repulsive color channels, that come with an 
opposite sign in the grand potential. It is tempting to conclude from this result that 
it provides an a posteriori justihcation of the success of approaches involving free 
quasiparticles in the description of the equation of state, even in the strongly coupled 
phase. 


An obvious drawback of our framework is the neglect of chiral symmetry, leading to results 
that may be inaccurate in the light quark sector. QCD in Coulomb gauge is currently the 
formalism which is maximally close to ours while fully including chiral symmetry. Some 
work has been done in the study of pure Yang-Mills theory and by using a toy model with 
conhning potential that mimics QCD [34]. Modelling the full quark-gluon plasma within 
Coulomb gauge QCD is however a huge task that still remains to be achieved. 

Some comments can be made about the large-iV behavior of our results. The meson 
masses depend on the number of color through the factor n,.qq = —4(1 — only, so 
the meson masses are of order 1 at large iV, with corrections in 1/iV^ as expected from a 
quenched potential. Moreover, it has been shown in [3] that the quark contribution to the 
equation of state behaves as Nf N in ’t Hooft’s limit, as expected. This is an important 
check of the ability of the present model to deal with the large-iV limit. 

Future developments of the present model should include the computation of the viscosity- 
over-entropy ratio. Such a computation can in principle be done without extra parameter. 
Hence, it is an important extension of our formalism that we hope to present in forthcoming 
works. 


33 



ACKNOWLEDGMENTS 


G. L. thank F.R.S-FNRS for financial support, and D. Cabrera, R. Rapp and C. Ratti 
for their interesting discussions and suggestions. 


[1] K. Yagi, T. Hatsuda and Y. Miake, Quark-Gluon Plasma: From Big Bang to Little Bang 
(Cambridge Monographs on Particle Physics, Nuclear Physics and Cosmology, Cambridge 
University Press, 2008). 

[2] D. Cabrera and R. Rapp, Eur. Phys. J. A 31, 858 (2007); Phys. Rev. D 76, 114506 (2007). 

[3] G. Lacroix, C. Semay, D. Cabrera, and F. Buisseret, Phys. Rev. D 87, 054025 (2013). 

[4] R. Dashen, S.-K. Ma, and H. J. Bernstein, Phys. Rev. 187, 345 (1969). 

[5] M. I. Gorenstein and S. N. Yang, Phys. Rev. D 52, 5206 (1995); A. Peshier, B. Kampfer, 
O. P. Pavlenko, and G. Soff, Phys. Rev. D 54, 2399 (1996); P. Leva! and U. Heinz, Phys. Rev. 
C 57, 1879 (1998). 

[6] F. Buisseret and G. Lacroix, Eur. Phys. J. C 70, 1051 (2010). 

[7] E. V. Shuryak and 1. Zahed, Phys. Rev. D 70, 054507 (2004). 

[8] G. Lacroix, C. Semay, and F. Buisseret, arXiv: 1408.0958 (2014). 

[9] J. 1. Kapusta and C. Gale, Finite-Temperature Field Theory. Prineiples and Applications, 
Cambridge Monographs (2006). 

[10] C. J. Joachain, Quantum Collision Theory (North-Holland, New-York, 1975). 

[11] O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, Phys. Lett. B 543, 41 (2002). 

[12] H. Satz, hep-ph/0602245; V. V. Dixit, Mod. Phys. Lett. A 5, 227 (1990). 

[13] O. Kaczmarek and F. Zantow, Phys Rev D 71, 114510 (2005). 

[14] G. S. Bali, Phys. Rep. 343, 1 (2001); C. Semay, Eur. Phys. J. A 22, 353 (2004); M. Cardoso 
and P. Bicudo, Phys. Rev. D 78, 074508 (2008). 

[15] S. Gupta, K. Hiibner, and O. Kaczmarek Phys. Rev. D 77, 034503 (2008). 

[16] A. Mocsy and P. Petreczky, Phys. Rev. D 73, 074007 (2006). 

[17] J.-P. Blaizot, E. lancu, and A. Rebhan, Phys. Lett. B 470, 181 (1999); Phys. Rev. D 63, 
065003 (2001). 

[18] S. Pratt and W. Bauer, Phys. Lett. B 329, 413 (1994). 


34 



[19] M. Jacob and G. C. Wick, Annals Phys. 7 , 404 (1959). 

[20] M. I. Haftel and F. Tabakin, Nucl. Phys. A 158, 1 (1970). 

[21] V. Mathieu, F. Buisseret, and C. Semay, Phys. Rev. D 77, 114022 (2008); F. Buisseret, V. 
Mathieu, and C. Semay, Phys. Rev. D 80, 074021 (2009). 

[22] G. S. Bali and A. Pineda, Phys. Rev. D 69, 094001 (2004). 

[23] J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

[24] O. Oliveira and P. Bicudo, J. Phys. G 38, 045003 (2011). 

[25] D. Binosi and J. Papavassiliou, Phys. Rep. 479, 1 (2009). 

[26] PDG, K. A. Olive et al, Chin. Phys. C 38, 090001 (2014). 

[27] P. O. Bowman et al, Nucl. Phys. Proc. Suppl. 119, 323 (2003); T. Hell, S. Rofiner, M. Cristo- 
foretti, and W. Weise, Phys. Rev. D 81,074034 (2010). 

[28] H. Satz, arXiv: 1501.03940 (2015). 

[29] Bazavov et al, PoS (Lattice 2013), 154; arXiv: 1312.5011 (2013). 

[30] S. Borsanyi et al, Phys. Lett. B 370, 99 (2014). 

[31] C. Ratti, private communication (2014). 

[32] J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34, 1353 (1975); E. V. Shuryak, Phys. Rep. 
61, 71 (1980). 

[33] S. Borsanyi, et al, JHEP 1208, 053 (2012). 

[34] P. M. Lo and E. S. Swanson, Phys. Rev. D 81, 034030 (2010); T. Yepez-Martinez, 
A. P. Szczepaniak and H. Reinhardt, Phys. Rev. D 86, 076010 (2012). 


35 



